Adiabatic quantum search with atoms in a cavity driven by lasers 
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We propose an implementation of the quantum search algorithm of a marked item in an unsorted 
list of N items by adiabatic passage in a cavity-laser-atom system. We use an ensemble of N identical 
three-level atoms trapped in a single-mode cavity and driven by two lasers. In each atom, the same 
level represents a database entry. One of the atoms is marked by having an energy gap between its 
two ground states. Appropriate time delays between the two laser pulses allow one to populate the 
marked state starting from an initial entangled state within a decoherence-free adiabatic subspace. 
The time to achieve such a process is shown to exhibit the Grover speedup yN. 

PACS numbers: 03.67.Lx, 32.80.Qk, 42.50.-p 



One typical problem of quantum computation concerns 
the search of a marked entry in an unsorted database by 
accessing it a minimum number of times. The Grover 
algorithm |l[ achieves this task quadratically faster than 
any classical algorithm. It is formulated in terms of a se- 
ries of quantum gates applied to a quantum register con- 
sisting of a collection of qubits encoding the database en- 
tries. An initial uniform superposition \w), independent 
of the searched state |m), is rotated step by step under 
the action of appropriate gates. The searched state is ex- 
hibited by an oracle function which checks if a proposed 
input is the searched state, returning for instance 1 in this 
case and otherwise. The number of steps grows as N 1 / 2 
with N the database size, whereas a classical algorithm 
requires on average N/2 calls. This quantum circuit al- 
gorithm has been tested experimentally for two qubits 
(N=4) by several techniques resting on NMR 0, (3[ , op- 
tics [3, Q and trapped ions Q. There have also been 
proposals of experimental implementations using cavity 
QED where the quantumgate dynamics is provided by a 
cavity-assisted collision [7( or by a strong resonant clas- 
sical field 0. 

A time continuous version of the Grover algorithm has 
been proposed by Fahri and Gutmann Q who, instead 
of using an explicit oracle, mark the searched state with 
an energy E while the others are degenerate with en- 
ergy 0, and use a driving Hamiltonian that leads con- 
tinuously the initial state to the marked one. Choosing 
an Hamiltonian V — E\w)(w\ to drive the free system 
Hq = E\m)(m\, they have shown that a Rabi-like half- 
cycle leads to the target marked state in a time growing 
as N 1 / 2 /E. Note that the Grover speedup is quadratic, 
independently of any increase of E with N which would 
simply amount to renormalizing the time. An experimen- 
tal realization of this analog Grover algorithm has been 
performed by NMR 10] in a setting where a quadrupo- 
lar coupling makes a spin 3 jl nucleus a two-qubit system 
(AT=4). 

Adiabatic versions of the time continuous Grover al- 
gorithm have been proposed 11, 13, IH mainly to take 



FIG. 1: Linkage pattern for the individual atoms. The un- 
marked atoms have two degenerate ground states \g) and \g ). 
One atom is marked with the state \g') shifted. The laser of 
Rabi frequency fi' (resp. f2) is resonant with the g' — e tran- 
sition for the marked (resp. unmarked) atom(s). The cavity 
of Rabi frequency G is resonant with the g — e transition. 



advantage of the robustness of adiabatic passage with 
respect to fluctuations of the external control fields as 
well as to the imperfect knowledge of the model. They 
have been formulated with a Hamiltonian which con- 
nects adiabatically the initial ground superposition |to) 
to the marked state \m) through an avoided crossing: 
H = [1 - u(t)]Hi + u(t)H f where H, = I - \w)(w\, 
Hf = I — \m){m\, and u(t) is a function of time growing 
from to 1. Roland and Cerf [l3| have shown that only 
a specific speed of the dynamics controlled by u(t) allows 
one to achieve the transfer to the marked state in a time 
growing as N 1 / 2 . 



In this paper, we show an implementation of the adi- 
abatic Grover algorithm based on a physical system, 
which is in principle scalable. This is, to our knowl- 
edge, the first proposed implementation of the adiabatic 
Grover algorithm. It is formulated with an Hamiltonian 
H = o + V(t) where Ho is considered as an oracle, and 
given, while V(t) is slowly varying in time and such that 
there is an instantaneous eigenvector adiabatically con- 
necting the initial superposition \w) to the marked state 
\m). We use an ensemble of N identical three- level atoms 
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trapped in a single-mode cavity of coupling frequency 
G, and driven by two lasers of Rabi frequencies f2 and 
f2'. The atomic levels are in a A configuration with two 
ground states \g) and \g') coupled to the excited state |e) 
by, respectively, the cavity and the two lasers (see Fig. [1]). 
This can be realized in practice by considering Zceman 
states and laser and cavity fields of appropriate polariza- 
tions. The states \g') of the N atoms are considered as 
the database entries. The energy of the state \g') of the 
marked atom is shifted by an amount 8 with respect to 
that of the unmarked atoms, which is set to zero. The 
states \g) allow for the coupling of all the atoms through 
the exchange of a single photon with the cavity (see Fig. 
EK). The initial state we shall start from is the entan- 
gled state \w) = \g',0) = (l/v^V)Ef=i IflJ.O) featuring 
a collective single-photon atomic excitation. Note that 
the label of each atom is added as a subscript 1, • • • , N 
and chosen so that the marked atom has tag N. Such 
a state can be prepared for instance before the marking 
of the atom using the stimulated Raman adiabatic pas- 
sage (STIRAP) technique 14 1, exactly as shown in [HI 
to store single-photon quantum states. 

The process we introduce here allows one to drive adi- 
abatically the population of the entangled state \g',0), 
which corresponds to the superposition of both the 
marked state |m) = |<7jy,0) an d an unmarked collective 
state \g' u , 0) that we introduce below, to the single marked 
state \g'jy,0) (see Fig. [2b). This process will be referred 
to as an inverse fractional stimulated Raman adiabatic 
passage (if-STIRAP) since it is a time inversion of the 
so-called fractional STIRAP (f-STIRAP) which transfers 
the population from a single state to a superposition of 
state [161]. This is implemented by first switching on Q 
and fi' together and next switching off tt 1 before fL The 
time to achieve such a process will then be shown to grow 
as y/W in order to satisfy adiabaticity. 



The Hamiltonian describing the system of N atoms is 

N 

H Q = S\g N ){g N \+u;J2\ e j)( e jl i 1 ) 
i=i 

We consider a cavity mode of frequency to and coupling 
strength G together with two lasers of frequencies oj, lu — 5 
and pulse shapes f2(i), fl'(t) which do not grow with N. 
The resonant driving provided by the atoms-cavity-laser 
system is described by 

iv 



(a) |e,v,0) |ei,0) |e 2 ,0) 



V = uja>a + Ga^2 \ e j)(9\ 



i=i 



N 



ne 1 ^ + fTe^-^J Wj) ( e j I + h - c - ( 2 ) 

The full Hamiltonian H = Hq + V has a photonic block 
diagonal structure. Each block is labeled by the number 
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FIG. 2: (a) Coupling scheme in the cavity. The cavity G, laser 
Q, laser fl' Rabi frequency are featured by respectively thick, 
thin, and dashed arrows, (b) Equivalent scheme where the 
states \g[, 0} (resp. |ei, 0)), i = 1, N — 1, of frame (a) form the 
collective unmarked ground state 13^,0} (resp. excited state 
|e u ,0)). The effective cavity Rabi frequency to the collective 
unmarked excited state is \/N — 1G. 



k of photons in the cavity when the N atoms are in their 
ground state \g). The corresponding multipartite state 
\gi ■ ■ ■ gn) ® \k) is denoted \g, k). All the states which 
are connected to \g,k) span a subspace whose projection 
operator is Pk- As each block is decoupled under H from 
the other ones, H — J2k=o PkHPk, we shall focus on the 
block PiHPi associated with a single photon in the cav- 
ity and show that it allows us to implement an adiabatic 
Grover search algorithm. 

The multipartite state \g, 1) is connected by to 
exactly two families of states as illustrated in Fig. [5^. 
Upon absorption of the cavity photon, the excited state 
|e) of any of the atoms, say atom j, can be reached 
while the other atoms remain in their ground state 
\g); the corresponding multipartite state is |ej , 0) = 
\9i ' ' • 9j-i e j9j+i • • • ) <8 |0). The state \ej, 0) can also be 
reached by absorption of one laser photon when any of 
the atoms, say atom j, is in the ground state \g') while 
the other atoms are in the ground state \g): \g'j,0) 

.'/! • ■ ■ gj~Wj9j+i ■■■qn)® |o). 

In order to remove the oscillatory time de- 
pendence introduced by the laser we consider 
atomic states which are dressed by laser and cav- 
ity photons, and use the resonant transforma- 
tion R = e- m \g' N ,0){g' N A + e^g, l)(g, 1| + 
Ef=i(l^0)<^,0| + e-- t |e J ,0)(e J ,0|). 

As we shall see, the states which are relevant for the 
Grover search are the N — 1 states \g'j, 0), which are un- 
marked, and the state |<^y,0) which is marked. Among 
the unmarked atoms, none should play a privileged role. 
Hence we shall consider them collectively and label the 
corresponding state with a subscript u. We rewrite this 
Hamiltonian in a new basis which features the uniform 
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superposition of the unmarked ground states 

N-l 

l^o> = -==£l^,o>, 



(3) 
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and the uniform superposition of the excited states asso- 
ciated with the unmarked atoms 



N-l 

|e«>0) = -_=^| ei) 0). 



(4) 
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This is achieved with the unitary transformation 
w = |^,0)( 3 5v,0| + | ejv ,0)( ejv ,0| + \g,l)(9,M + 
EZ=\ U ld(Wi,0){9' j ,0\ + \ei,Q}{e j ,0\) where U is any 
unitary matrix such that the (N — 1) elements of one 
of its columns are equal. In this new basis, the part of 
the Hamiltonian restricted to the subspace spanned by 
the states |^,0), |^,0), \g, 1), |e u ,0), \e N ,0) is decoupled 
from the rest. It reads 



Hi = 



Vt 

v} v c 



(5) 



with 
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and £ = fi + e-^'ST, £' = n' + e^ft. This Hamiltonian, 
whose derivation is exact, is represented in Fig. (2b- 

By means of a unitary transformation T, we diagonal- 
ize the block V c which admits the eigenvalues 70 = and 
7± = ±\/~NG whose associated eigenvectors read 



•n) = -i|e w ,o)--^| e „,o) 



(6) 



l7±) = ^(^v |e - 0) 



1-I|e u ,0)±| 5 ,l) 



The new Hamiltonian reads, 
K,0),|^,0},| 7 o},l7+},l7-) : 



with 
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in the basis 
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The time evolution of the non-resonant components in 
£ and £' is much faster than the evolution of f2 and f2' 
which occurs over a time scale T S> <5 _1 . Hence it is 
justified to replace these contributions by their vanishing 
average values over times r'«T«T: £ ~ SI and £' ~ 
SI', where f(t) = i f* +r duf(u). This is the resonant 
approximation. Similarily, the unitary evolution of C 
is much faster than that of A if Sl pca k/-/VG <§; 1 since 
the respective eigenvalues scale as VNG and Sl pca k/-\//V 
where Sl poa k is the peak amplitude of the pulse SI. Upon 
performing an adiabatic elimination we thus obtain an 
effective Hamiltonian H e g — A~BC~ 1 B^ which contains 
all the contributions up to order (Sl pea k/-/VG) 4 . In the 
basis \g' N ,0), I70), |ff„,0), the effective Hamiltonian reads 
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Our aim is to transfer adiabatically the population 
from an initial state \g' , 0) which gives no privileged role 
to any of the N states \g'j, 0) to a final state which coin- 
cides with the marked state \g' N , 0) in a time which scales 
as V^V- The population transfer mechanism is most eas- 
ily revealed in the basis of the instantaneous eigenstates 
o{H eS (t) 

\0)(t) = cos8(t) \g' N ,O)-sm0(t) \g' u ,0) (10) 
|±>(t) = -^=(sm9(t) \g' N ,0)+co S 6(t) |^,0)±| 7o )), 

pertaining to the eigenvalues and ±A(£) where 

A(t) = ^=J(N-i)n> 2 (t) + W(t). (11) 



N 

The instantaneous angle 9{t) is defined through the rela- 
tion 



tan„( f ) = -^m- 



(12) 



Requiring the instantaneous eigenstate (fH))) to coincide 
at the initial time ti with the uniform superposition 



\g',0) = -l=\g' N ,0) + ^l-l-\g' u ,0), (13) 



and at the final time with the marked state \g' N , 0) entails 
that 



tan0(ti) = -\/N - 1, tan0(£ f ) = 0. (14) 

This implies that the two pulses must be switched on 
simultaneously, SI '(ti) — fi(ti) and that the pulse SI' is to 
be turned off before SI. In the adiabatic representation 
(fT0|) . the Hamiltonian (J9j) reads 



A 71 £ 







(15) 



V2 L 



-A 
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where = 1+t \ n i g ^ tan 6. In the adiabatic regime, the 
transitions between instantaneous eigenstates are negli- 
gible. This will be achieved if the Hamiltonian varies 
sufficiently slowly in time so as to keep <C A. On the 
other hand, we wish to control the process duration and, 
in particular, to prevent it from becoming arbitrary large. 
For that purpose, as proposed in [l3j], we choose to re- 
quire 9 and A to be in a constant (small) ratio e at all 
times, independently of N: 



= eA. 



(16) 



This choice is adopted for the sake of clarity, since in 
this case the scaling can be analytically determined and 
proven to scale as s/~N. We have investigated the robust- 
ness of this approach. 

Given a laser pulse f2, this equation will allow us 
to determine the pulse f2' which is needed to remain 
in the instantaneous eigenstate \0}(t) with a probabil- 
ity larger than 1 — e 2 throughout the process, starting 
from the uniform superposition |0)(ti) = \w), and end- 
ing up in the marked state |0)(if) = 1.9^,0) after some 
time T = tf — t{. Indeed, rewriting (fTTj) with (fl2|) as 
a = -Wi + tan 2 en, we obtain from (1161) a differential 

V N 

equation for tan#, i. e., for the ratio f2'/J7. Its solution 
satisfying the initial condition (fl4|) reads 



i - 



sA(t) 



(17) 



1 + eA(t){2y/W=T - eA{t)} 



where we define A{t) = J duQ(u). The process duration 
T is obtained implicitly upon specifying that at time tf 
the ratio on the left side of (fT?)) vanishes: 



eA(T) = y/N - 1. 



(18) 



Expressing the total area of the pulse fl in terms of its 
average amplitude £1, A(T) = f2T, we arrive at 



(19) 



This shows that, for a constant average amplitude (i.e. 
independent of N), the duration scales as y/~N \ Note 
that we can equivalently increase the average amplitude 
as ^/N for a constant time T. We can determine from 
(HTJ) that VtT grows as (</N-l)/y/N - 1, i.e., asO(N°). 




Time t (x 1/T) 

FIG. 3: Numerical dynamics with the effective Hamiltonian 
([9]) for N = 8, e = 0.05 and a Gaussian Rabi frequency 
n(t) = n pC ake- (t/T)2 with fipeakT = ^iV^T/(e0F) [in order 
to satisfy (TJHJ) and (TH)])]. fi'(t) is determined from |T7J). Top: 
Rabi frequencies. Bottom: Populations of the collective un- 
marked state P u (t) = \{g' u ,0\(j))(t)\ 2 , and of the marked state 
Piv(f) = \{g'ff, 0\4>)(t)\ , where \4>){t) is the dynamical state 
vector. 



|0)(t) in the adiabatic limit. Notice that the choice (jT6j) , 
which leads to the seemingly complicated pulse relation 
(fT?)) . gives in practice a simple smooth bell-shaped pulse 
(see Fig. [3]). We have numerically checked that the effi- 
ciency of the transfer is, as expected, preserved for a wide 
range of N with f2 pea kT growing as y/N and an almost 
constant n' pcak T. 

In conclusion, we have proposed the first physical 
implementation of the adiabatic Grover search using a 
cavity-laser-atom system and robust processes related to 
STIRAP. The calculation has been conducted with pulses 
based on the constraint (fT6]) that has allowed us to prove 
the scaling analytically. We have checked the robustness 
of the \/N scaling by numerical simulations using other 
less restrictive adiabatic pulse shapes satisfying (fT4| . 

The authors are grateful to N. J. Cerf and H.-R. Jauslin 
for useful discussions and acknowledge the support from 
the EU projects COVAQIAL and QAP, from the Bel- 
gian government programme IUAP under grant V-18, 
and from the Conseil Regional de Bourgogne. 



Figure [3] displays the pulses and the population dy- 
namics resulting from (fT7|) with N — 8, e = 0.05 and a 
Gaussian pulse Q of characteristic duration T. As pre- 
dicted, the transfer to the marked state is very efficient. 
The very low transient population in the excited states 
stems from the fact that the dynamics is expected to 
remain in the instantaneous decoherence-free eigenstate 
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